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Abstract 

We use a scheme of separation of degrees of freedom for a system in order to 
produce two systems with finite number of degrees of freedom. Our intent is 
to measure the energy square relative fluctuation (SRF) of the observable part 
through the simulation of two simple examples of composed systems, the simple 
harmonic oscillator, and the chain of quartic oscillators. We want to test the 
result found previously by us through the finite heat bath canonical ensemble 
(cond-mat/0210525), which is an application of Tsallis' statistics. We see that 
the results obtained here are in very good agreement with the theoretical pre- 
dicted values. This suggests that this kind of finite systems are ergodic, and 
that they do not provide "bad" statistics. 
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1 Introduction 



Tsallis' statistics, first proposed 14 years ago J], is a field which attracts great 
attention nowadays. It has been applied to a variety of fields, among which 
are included: fully developed turbulence 0, anomalous diffusion jSj, and the 
Nose- Hoover thermostat 

Certainly, many proposals to explain its physical meaning, in particular the 
interpretation of the generalization parameter q, have risen. Some of them 
|S] suggest that q measures the fluctuation of a quantity pertaining the phc- 
nomenom in which the interpretation was proposed. 

We follow the path initially proposed by Plastino and Plastino 6 , and Almeida 
[7]. This approach is based on the idea that Tsallis' power law canonical distri- 
bution is the probability function of a system which is in thermal equilibrium 
with a heat bath which has a finite number of degrees of freedom, and the quan- 
tity is the exponent of the energy in the accessible phase space volume of 
the z-th system which is proportional to the number of degrees of freedom of 
this system. It may either be positive or negative. 

It should be noticed that in this formalism the distribution is obtained in the 
first place (instead of postulating an entropy functional and then finding the 
canonical form by optimization of the former) from the system's phase space 
geometry. 

Of course, this is only one case where Tsallis' statistics can be said to "apply" , 
since the generalization parameter q, firstly free of interpretation, provides a 
very convenient way of describing a wide class of phenomena. Henceforth, we 
will refer to our framework as the finite bath canonical ensemble, but still use 
the notation introduced previously. Any quantities labeled by 0, 1, and 2 belong 
to the composite (observable+bath) , observable, and heat bath system. 
Many theoretical developments were achieved along this line [HlElEmj- Also, 
a few simulations were performed in order to test the present approach |12[I13| . 
A detailed description of the finite heat bath canonical ensemble may be found 
within these references. 

We propose here to simulate two simple systems in order to measure the square 
relative fiuctuation (SRF) of a system. Wc want to test the result found in jl4| . 
We do this by using a procedure in which we separate the degrees of freedom 
of a classical system, which is not is thermal contact with any heat bath, and 
define one of them as the observable system, and the other as the bath. This 
was primarly used in |12[ I13j . 

We organize the paper as follows: in section 2, we give a brief explanation of the 
relative fiuctuation calculation made in ;14[. Section 3 is devoted to introduce 
the model, the systems, and the simulations' results. In section 4 we discuss the 
physical meaning of this approach. In section 5 we address our conclusions. 
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2 The relative fluctuation in the finite bath ap- 
proach 

In order to calculate the SRF of a random quantity x, we need to calculate the 
following expression: 

The angular brackets denote ensemble averages. 

In J2] , the authors calculated equation for a system's energy Ei . This energy 
is distributed according to the following power law canonical distribution: 

= (2) 
^1 

where Zi is the partition function, the parameter (3* satisfies: 

rfe-l) = ^, (3) 

El and Eq — Ei E2 are the observable and total energies. Equation Q is 
just another way of writing the distribtuion function of a system in contact with 
a finite heat bath, therefore, it is not, from the very beginning, a distribution 
which arises due to the presence of unusual characteristics, such as long-range 
interactions. All the moments of Ei were calculated in ^3], and equation |^ 
was written as: 

{AEl) _ q,/{q, - 1) 1 



(4) 



{Elf - 1) go/(9o - 1) + 1 

where: 

^ = ^ + (5) 
90-1 (71-1 g2 - 1 

This expression is nothing less than the additivity of the exponents of the al- 
lowed phase space volumes of the three systems. 

The temperature of a system was defined in |9j , through the equipartition the- 
orem, and is given by: 

k T = ^0 = (^1) ^ {E2) ... 
90/(90-1) 91/(91-1) 92/(92-1) ^' 

It may seem confusing the 3 q's introduced above, one for each system, since 
the whole non-extensive statistics was primarily based on a single parameter q. 
This approach necessarily leads us to the 3 g's, since each system has its own 
phase space, and therefore, characterized by its own q. Fortunately, only one is 
really important, namely 52- This is the number which characterizes the distri- 
bution function of the observable system, and if this q goes to 1, inother words, 
if the number of degrees of freedom of the bath goes to infinity, the distribution 
((2)) goes to the Boltzmann-gibbs (BG) exponential, and the whole formalism. 
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to the classical BG one. 

The simulation performed here intend to test the validity and accuracy of 
The results are presented in the next section. 



3 Numerical models 

The idea here is as follows: suppose we have a system of N particles in D 
dimensions isolated from any other kind of interaction, which means that it is 
not coupled to an external bath. It has the following hamiltonian: 

DN 2 

H = J2^ + U{x,) + U,nt{x,), (7) 

2—1 ^ 

where the first term in the right is the kinetic energy. The other two represent 
the potential part. Note that there may be an interaction term in this system. 
Then we pick one degree of freedom, say pi , and separate it from the rest of the 
hamiltonian: 

2 DN 2 

= TT^ + V + Uix,) + Um{x,) =Ei+ E2. 

^ i=2 * 

Now we consider pi as our observable system and the rest of the original hamil- 
tonian as the heat bath. Since the number of degrees of freedom of Ei is finite, 
in this case equals to 2DN — 1, provided that DN is finite, the distribution of 
El , or equivalently pi , is always a power law, and the results found in references 
IHlElIinillllinilSlini can be tested. 

Although it seems to be very simple to work in this approach, it should be 
employed very carefully. Some cases are not direct as they seem to be, we must 
take into account the mixing of the coordinates. 

A great example is provided by the three dimensional harmonic oscillator, where 
Uint{xi) = 0. In this system we have three momenta and three positions. Hence, 
it is tempting to use equation to pi considering all the other degrees of free- 
dom as the heat bath. In this case we would have the following numbers: 

go ^ g. gi ^ 1 . g2 ^ 5 
go-i 'gi-l 2' 92-1 2' 

These are not correct by the fact that pi only trades energy with its correspon- 
dent position xi, the other coordinates (p2 : Pa j 2:2 , ) do not participate as a 
bath. In this case, the correct quantities are: 

go gi _ 1. g2 _ 1 /gN 

go-1 'gi-1 2'(?2-l 2' ^' 

For this reason, we cannot consider the sum of any momenta as a system, for it 
will not span the whole allowed phase space. Therefore, it is not ergodic. 
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In a system which Uint{xi) ^ 0, all coordinates trade energy among themselves, 
and the bath for a single, or more, degree of freedom is defined by the rest of 
the coordinates. Our second simulation is one of this case. When the system is 
coupled to an external bath, this carefuUness is not needed since this interaction 
makes all the degrees of freedom to span the whole allowed phase space. 
This approach was inspired by Khinchin ^S] who suggested that degrees of free- 
dom could be separated as a system and bath. 

Thus, the simulations we present here do not describe a physical particle. They 
describes only a degree of freedom of a system. 



3.1 The single harmonic oscillator in one dimension 

The first system to be simulated is, perhaps, the simplest one to consider. Its 
hamiltonian reads: 

H=^ + -Kx^. (9) 

We take Ei as the kinetic part, and the bath as the potential term. The values of 
q for this problem are given in (jSJ. We notice that this is a symmetric situation, 
in the sense that if we choose the potential part as the observable system, the 
results will be the same. The SRF for this system reads: 

— y- = 77- (10) 

The result of the simulations is presented in figure 1. It should be noticed that 
the SRF does not deviate much from the predicted value. 
As a by product of the simulation, we present, in figure 3, the fit for the 



4e+05 6e+05 
Simulation time (t) 



Figure 1: Relative fluctuation of the kinetic part of the one dimensional oscillator 
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distribution of Ei. The slope values —0.481. Theoretically predicted value is 
given by the power of the structure function of the momentum: 

= -0.5. (11) 

In this simulation, we have taken the following values for the important pa- 

6 I ' , ' , ■ , ■ 1 




2 I , 1 , 1 . 1 . 1 

-6 -4-2 2 

ln(1-E,/E„) 

Figure 2: Linear fit of the distribution of Ei 

rameters: mass: m = 0.5, elastic constant: K ~ 2.0, time step: At — 0.01. The 
data for building the fit of the distribution were taken after the thermalization 
time was achieved. We estimated it at about: tther = 7.0 x lO''. We did this by 
observing the time evolution of Q{t) — Xi^j--^, where Xi represents the canon- 
ical variables, Q is proportional to the ratio j , Vq{Eq) is the composed 
system's accessible volume and VloiEo) is its structure function, and is given by: 



qo/iqo - 1)' 

and the overbar denotes time average. We test the approach of it to unity, for 
this can be a good estimate of the thermalization time of the system. 
We used a Leapfrog routine to integrate the system. 



3.2 The chain of quartic oscillators 

This system was inspired by the Fermi-Pasta-Ulam model, and was first studied 
in ^21- Its hamiltonian is given by: 

(12) 




i=l 
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where {pi, Xi) denote the momentum and position of an oscillator, and N is the 
total number of oscillators. 

First we will measure the SRF of the kinetic energy of the first oscillator. We 
choose 
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N 2 ^ 
i=2 i=l 

From a simple scaling argument jl2| . we can easily calculate the g's for this 
problem. They are given next: 



90 



3iV. qi 1. 92 

4 2' 92-1 



3iV 



9o - 1 

Therefore, we have a relative fluctuation given by: 



3A^ + 4 



(13) 



In figure 4, we present the results for several values of N. 
In table 1, we show the SRF for the values of N considered in the simulations. 
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5e+06 1e+07 

Simulation time (t) 

Figure 3: Relative fiuctuation of (pi) for iV = 8, 12, 16, 32, 64, 128. The lowest 
line represents the first case, and the others follow the sequence. 

We see that while the number of oscillators increases, the SRF approaches 2, 
its iV — > oo value. 

We also considered a system composed of momenta pi,p2, only, and then one 
composed of the first two plus ps, and observed that the SRF agree with (|13(l . 
This is to be expected since the canonical coordinates are all connected by the 
quartic interaction. 
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N 




08 


1.571 


12 


1.7 


16 


1.769 


32 


1.88 


64 


1.938 


128 


1.969 



Table 1: Relative fluctuation of ^p^) for several values of N 



As a second test, we studied the average value of the potential part of this 
hamiltonian, namely: 



N 



i=l 



{Xi+i - Xi)'^ 



(14) 



The number of degrees of freedom for this system can be directly calculated by 
observing that for the kinetic part (which plays the bath's role now), we have: 

92 N 



92-1 2 ' 

then, we have for the potential energy (observable system): 

qi ^ N 
91 - 1 4 ■ 

The SRF, equation I^J is given by: 

{AEif 8 



(15) 



(16) 



In figure 4, we present the time evolution, after thermalization, of (|16() for several 
values of N. In table 2, we calculate the values of the relative fluctuation for 
each case. Once again, if the number of oscillators is increased, the SRF goes 
to the N oo value, i. e., zero. 

A third test of relation © would be to study only a term in the potential 
energy, say Xq. The 91/(91 — 1) values may be, in some cases, not so hard to 
calculate, and Q would be readly applicable. Unfortunately, this procedure 
does not work, as we observed by studying the above term. Although its mean 
value and SRF converge, these values are not predicted by our calculations. We 
observed that this term does not obeys the equipartition, in other words, its 
mean value does not follow: ^ 

(4) = -^keT. 
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Figure 4: Relative fluctuation of the potential energy for A'' = 8, 12, 16. The 
highest line represents the first case, and the others follow the sequence. 



N 


(AE-i) 


08 


0.285 


12 


0.2 


16 


0.153 


32 


0.08 


64 


0.04 


128 


0.02 



Table 2: Relative fluctuation of the mean potential energy for several values of 
N 



This suggests that isolated terms in the potential are non-ergodic, although the 

whole potential energy is an ergodic one. 

In figure 5, we show the time evolution of the quantity: 

which, by the former, should approach 0.25. Clearly, we see that it does not. 
We have also observed this behaviour in other terms, such as ^\ ' 
Although the cquipartitioii docs not hold in this case, we could use the final value 
of 91/(^1 — 1) and calculate the relative fluctuation of this term. Once again, the 
value obtained does not coincide with the one observed in the simulations, not 
shown here. Hence, we cannot predict, with the tools of the canonical ensemble, 
for this Xq term, any quantity, and, by the nature of the problem, we can also 
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say that the other terms in the potential energy follow this line. 
The reason for this discrepancy may lie in the fact that no term in the potential 
energy is independent of any other term in the potential energy. In the kinetic 
part, no single momentum component is free, it interacts with the potential 
part, but it does not interact with the other momenta components. This is why 
the single terms in the potential energy do not span the allowed phase space. 

0.12 I ^ , ^ , ^ , ^ , ^ 1 



0.115 - 




0.095 

0.09 I ' ' ' ' ' < ' ' ' 

2e+Q6 4e+06 6e+06 8e+06 1e+07 

Simulation time: t 

Figure 5: Time evolution of iov N = 8. 



4 Discussion 

We saw that our separation of degrees of freedom scheme provides a simple way 
of studying systems with a finite number of degrees of freedom. The results 
obtained here are in very good agreement with the theoretical results. This 
fact indicates that the time averages of such finite systems coincide with their 
ensemble averages. This seems to be in contradiction with the thought that 
for statistical mechanics to apply, it is required infinite number of particles, or 
infinite number of degrees of freedom. We state that this represents no contra- 
diction at all for ensemble averages are taken with respect to the system's states 
of energy, and not with respect to the system's degrees of freedom. The num- 
ber of states of energy of a system are always infinite in number for a classical 
system (which is the case here). 

For example, consider a classical particle in a box. How many states of energy, 
i.e. how many values of the canonical coordinates {p,x), exist such that are 
compatible with the particle's energy? Clearly infinite. Generally, every degree 
of freedom defines a surface in its own phase space, with an infinite number of 
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points (states). 



5 Conclusions 

We have simulated two simple systems not in thermal equilibrium with a heat 
bath, in order to separate its degrees of freedom, and consider them as separate 
systems, one as the observable and the other as the heat bath. We intended 
to study systems with finite number of degrees of freedom, and therefore, to 
test the results found in ^1] for the square relative fluctuation of the observable 
system energy. 

Our results are in very good agreement with the theoretical predicted value, 
hence indicating that these systems are indeed ergodic. 

We can conclude that these finite systems do not represent a case of "bad" 
statistics, in the sense that the methods of statistical mechanics do not apply 
to them. Certainly, the method of Boltzmann-Gibbs, which assumes an ideal 
thermostat in certainly not appplicable, but our procedure of considering the 
system decoupled from the ideal heat bath, hence not in equilibrium with such 
a system, is clearly a good framework to use. 
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